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Abstract. We study initial-boundary value problems for the La- 
grangian averaged alpha models for the equations of motion for 
the corotational Maxwell and inviscid fluids in 2D and 3D. We 
show existence of (global in time) dissipative solutions to these 
problems. We also discuss the idea of dissipative solution in an 
abstract Hilbert space framework. 



1. Introduction 

Recently, so-called alpha-models of fluid mechanics have attracted 
attention of many researchers, see e.g. [3l|8l|9l[l5l[T6l[l9l|20l[2Tl|22l 
2M EH] and references therein. These models turned out to be rather 
ubiquitous, being relevant in such issues as turbulence and large eddy 
simulations, and, on the other hand, being related to the second grade 
fluids. 

In this work, we study the initial boundary value problem for the 
corotational Maxwell-alpha viscoelastic fluid flow: 

(1-1) + 2^"^^ + >^^.Vw. + Vp = Diva, 

i=l * i=l 

[1.2) a + X\^ + y ]ui^ + aW -Wa] =2r]S, 



;i.3) 

:i.4) 
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(1.5) U =0, 

an 

(1.6) u\t=Q = a, a\t=o = ctq- 

Here, is a domain in the space M"", n = 2, 3, m is an unknown 
velocity vector, p is an unknown modified pressure, a is an unknown 
deviatoric stress tensor, v is an auxiliary variable (all of them depend 
on points x in the domain Q, and on time t); 

£ = £{u) = {£,,{u)l £,^{u) = U^ + ^''^ 



W = W{u) = {W,,{u)), W.,,{u) 



2 \dxj dxi 

is the strain velocity tensor, 

1 / dui duj 

2 \dxj dxi 

is the vorticity tensor, 77 > is the Maxwellian viscosity of the fluid, 
A > is the relaxation time, a > is a scalar parameter, a and (Tq 
are given functions. The external force is, for simplicity, assumed to 
be zero. 

Equation (11. ip is the Lagrangian averaged (Euler-alpha-like) Cauchy's 
equation of motion, (11. 2p is the corotational Maxwell constitutive law 
P3j . (II. 4p is the equation of continuity, and (II. 5p is the no-slip condi- 
tion. When a = 0, (II. ip - (II. 6p becomes the initial-boundary value 
problem for the equations of motion for the corotational Maxwell vis- 
coelastic fluid [HI [231 [3ll SH] . Thus, (JTI]) - ([L6]) can be considered as 
an appropriate a-model for the corotational Maxwell fluid. 

In the particular case 17 = and (Jq = 0, we recover the Dirichlet 
problem for the celebrated Euler-a model: 



(1.7) 

(1.8) 
(1.9) 

(1.10) 
(1.11) 



o n r~i n 

OV \ ^ OV \ ^ „ „ „ 



dt ^ dxi 

i=l 1=1 



V = u — a'^Au, 
divu = 0, 



u 



0, 



an 



u\t=o = a. 



The Euler-alpha (also known as Lagrangian averaged Euler and invis- 
cid Camassa-Holm) equations were introduced and derived in [20l [2T] . 
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They are well-posed on small time intervals |3H1 E?]. However, nei- 
ther strong nor weak general global solvability result has been known 
for the three-dimensional domains. A sort of a Beale-Kato-Majda cri- 
terion for these equations was proposed in [22]. Some reviews of math- 
ematical results on the Euler-a model and many references may be 
found in [3l[8l[T5l|2l]. 

The Maxwell model is one of the basic and classical models of a vis- 
coelastic material. Its mechanical analogy is comprised of a spring and 
a dashpot connected in series [31]. The multidimensional Maxwell mod- 
els generate complicated systems of PDEs due the frame- indifference re- 
strictions and consecutive involvement of objective derivatives [2^ 
The simplest objective derivative is the corotational (Jaumann) one, 
in the incompressible case its use yields f ll.ip -f fT74|) with a = 0. 

Very few mathematical results are known for the corotational Maxwell 
fluid equations (see [HI |3ll |33]). In particular, there is no global solv- 
ability theorem, even in 2D. Moreover, there are evidences of non- 
existence of smooth solutions [5^ [T5]. Let us also mention paper [5B] 
with some numerical issues regarding the corotational Maxwell fluid. 

In these circumstances, if we want the problem to be solvable glob- 
ally, a possible way out is to consider a kind of a generalized solution 
different from the standard hydrodynamical weak solution framework. 
We are going to use the concept of dissipative solution due to P.-L. 
Lions. It was suggested in [26] for the Euler equations of ideal fluid 
flow, which have only been proven to be globally weakly solvable on 
the torus (this is a very recent result [16]). Some properties of dissipa- 
tive solutions are discussed in [HES]- Later, existence of such solutions 
was established for Boltzmann's equation [2S] (see also [13); and for 
various models arising in magnetohydrodynamics [ITJ [2], diffusion in 
polymers [15], and image restoration |40] . 

The notion of dissipative solution plays a key role in the problem 
of transition from kinetic theory to hydrodynamics. In the Euler hy- 
drodynamic limit, the renormalized solutions of Boltzmann's equation 
tend to dissipative solutions of Euler's one [Sj [351 112]. Other issues con- 
cerning relevance of dissipative solutions for the Euler equations and 
relation of this concept to the weak and measure- valued solutions may 
be encountered in [B] HO] ■ 

In math literature, the expression "dissipative solution" has various 
meanings. In particular, the notion that we use differs from the ones 
from [71[II1 [2811291 ED]. 

The objective of our paper is to introduce dissipative solutions for 
the corotational Maxwell-alpha problem (11. ip - (II. 6p . and to show their 
existence and basic properties. These solutions are always global in 
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time. In the appendix to the paper, we present the skeleton of the idea 
of dissipative solution via considering it in an abstract Hilbert space 
setting. 

Throughout the paper, for definiteness, we assume rj > 0. However, 
the results remain valid in the Euler-alpha case, and the proofs are 
similar but simplified (see also Remark 12. ip . 

Our approach is not directly applicable to the corotational Maxwell 
fluid (a = 0), but it allows us to construct some "ultra-generalized" 
solutions for that model (see Remark 13. 3p . 

The paper is organized in the following way. The next section con- 
tains preliminary material which is required for the formulation of the 
main result (Theorem l2.ip . The proof of the theorem is provided in the 
third section, which also contains some discussion of related open prob- 
lems. The appendix with a general explanatory approach to dissipative 
solutions is rather an illustration than a collection of results, therefore 
its framework is not general enough to encompass the Maxwell-alpha 
or Euler equations. 

2. Preliminaries, notation and the main result 

By M*^^", we denote the space of n x ra- matrices with the following 
scalar product: for A = {Aij), B = (Bij), 

n 

A : B = ^ ^ AijBij, 

and let M^^" be its subspace of symmetric matrices. 

Below in the paper, Q is considered to be a domain (i.e. an open set 
in M", n = 2, 3) possessing the cone property. We recall [1] that this 
means that each point x G is a vertex of a finite cone contained 
in Q, and all these cones are congruent. A finite cone is a set of the 
form 

= B,n {x + ay - x)\y e B2,^ > 0} 

where Bi and B2 are open balls in R", Bi is centered at x, and B2 does 
not contain x. We require this property to have Sobolev embeddings. 
We use the standard notations Lp{n,F), W^{n,F), H^{n,F) = 

W^in, F), H^{n, F) =W 2 (^, F) (/3 > 0) for the Lebesgue and Sobolev 
spaces of functions with values in a finite-dimensional space F. In this 
notation, sometimes we only keep the function space symbol and omit 
n and F, especially when F = Mg''". 

The Euclidean norm in finite-dimensional spaces F is denoted as | ■ | . 
The symbol || ■ || will stand for the Euclidean norm in L2{^l,F), and 
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II ■ 11/3 will stand for the same thing in H^{Q,F). The corresponding 
scalar products are denoted by (■, ■) and (■, O/S? respectively. 

Let V be the set of smooth, divergence- free, compactly supported 
in Q functions with values in M". The symbols H and V denote the 
closures of V in L2 and H^, respectively. We also use the spaces Vi = 
V n H\ i = 2,3, with the scalar product inherited from H^. 

There exists a continuous Leray projection 

P : L2in,M.^) H. 

Let us introduce the operator 

Aa = I - a^A, 

where a is the same as in fll.3p . and / is the identity map. 

In the space V, along with the scalar product (■, Oi^ "^^ use another 
one 

{u,v)v = {u,v) + (aVu,aVf), 

and the corresponding Euclidean norm || ■ ||y. Let us mention the 
inequality 

||u||i < max{l, l/a;}||M||v. 

We recall the following abstract observation [HI US] . Assume that we 
have two Hilbert spaces, X C Y, with continuous embedding operator 
i : X ^ Y, and i{X) is dense in Y. The adjoint operator i* : Y* ^ X* 
is continuous and, since i{X) is dense in Y, one-to-one. Since i is one- 
to-one, i*{Y*) is dense in X*, and one may identify Y* with a dense 
subspace of X*. Due to the Riesz representation theorem, one may 
also identify Y with Y*. We arrive at the chain of inclusions: 

(2.1) XCY = Y*CX*. 

Both embeddings here are dense and continuous. Observe that in this 
situation, for / G F, u G X, their scalar product in Y coincides with 
the value of the functional / from X* on the element w G X: 

(2.2) (f,u)y = {f,u). 

Such triples (X, y, X*) are called Lions triples. 

We use the Lions triples (\/3, V, 1^3*) and {H^, L2, H]^"^). In the first 
triple, the structure on V is determined by the scalar product {■,-)v- 
In the second one, we write H]^^ for {H"^)* due to relation of this space 
to Neumann boundary value problems |44j . 

The symbols C{J'; E), Cw{J] E), L2{J; E) etc. denote the spaces of 
continuous, weakly continuous, quadratically integrable etc. functions 
on an interval C M with values in a Banach space E. We recall that 
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a function u : J ^ E is weakly continuous if for any linear continuous 
functional g on E the function g{u{-)) : J' M. is continuous. 
We require the following spaces 

Wi = W,{Q,T) = {re L2{0,T;V3), r' G ^2(0, T; V^;)}, 

\\t\\wi = ||T|U2(0,T;y3) + lk'IU2(0,r;V3*), 

W2 = W2{n,T) 

= {re L2(0, T; H\Q, M^x")), t' G ^^(O, T; H;,\Q, RTI)}, 

\\'t\\w2 = lklU2(0,r;H2) + lk'llL2(0,T;H-2)- 

Here and below the prime symbol stands for the time derivative. Note 
that 

W,CC{[0,T];V), W2CC{[0,T]-L2), 

(2.3) ^(u,u)v = 2{u',u), ueW,, 
at 

(2.4) j^{u,u) = 2{u',u), ueW2, 
and 

(2.5) ^(u, 0)v = {u', <p), ueWu<peVs. 
at 

(2.6) ^(w,0) = «0), ueW2A^H', 

these facts are consequences of e.g. [HI Lemmas 2.2.7 and 2.2.8] and 
formula (12. 2p above. 

Let us introduce the operators 

and 

A2:H^^H^\ (^20-,$) = ((T,<l>)2, 
where (p and $ are arbitrary elements of V3 and H^, respectively. 

Let fi = rj/X. Consider the following formal expressions, where w 
and T are vector- and matrix-valued functions of time, respectively, 
and 5 is a positive number: 

Ei{w,t) 

(2.7) = -P^2^^^ - pj^i^/^M^^w, + PDivr, 

1=1 * i=l 

E2{w,t) 
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(2.8) = -T - 7^ - 5Z ^^TH ~ ^^^""^ + ^^""^^ + 2^^^"^^ 



A dt ^ dxi 

Ei{w,T,S) 

(2.9) = - 5P X] w,^^ - 6P f^i^a^h"^^'^ + Div r, 

i=l * i=l 

E2{w,r,6) 

(2.10) = - ^ - ^Yl ^^IT' ~ -^^^H + ^W{w)t + 2(5/i^(u;). 

i=l 

The symbol C will stand for a generic positive constant that can 
take different values in different lines. We shall sometimes write Cq to 
specify that the constant depends on the domain Q only. 

Let us recall the well-known inequalities 

(2.11) \\uv\\ < Cn\\u\\2\\v\\, u e H'^,v G L2, 

(2.12) II < Cf^||M||i||t;||i, u,v G 

(see e.g. |1H1 Corollary 2.1.1]). 
The following Gronwall-like lemma will be useful. 

Lemma 2.1. Lemma 3.1]^ Let f,x,L,M : [0,T] R be scalar 
functions, XjL,M G Li(0,T), and f G Wi{0,T) (i.e. f is absolutely 
continuous). If 

(2.13) x{t)>0,L{t)>0 
and 

(2.14) f'{t) + xit)<L{t)fit) + Mit) 
for almost allt G (0,T), then 

t 

fit)+ [ xis)ds< 



(2.15) exp L{s)ds 

for allte [0,T]. 
Now we can give 





/(0)+ /exp I / L{i)di I M{s)ds 
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Definition 2.1. Let a e V, ao e L2(ri, M^"""). A pair of functions 
(m, 0") from the class 

(2.16) ueC^{[0,ooy,V), aGC^([0,oo);L2(n,M^^")), 

is called a dissipative solution to problem (11.11) - (11.61) if, for all test 
functions C G C\[0, oo); V3), 6 G ^^([0, 00); i^2(fi,M^^")) and all non- 
negative moments of time t, one has 

2i,\\u{t) - ami + Mt) - 9{t)r 

<exp I / T{s)ds'^ |2;,||a-C(0)||^+||ao-^^(0)f 



(2.17) +2{E,{C,9)is),ais)-9is))]ds 
where 

r(t) = 7max{l, l/«^} 

2n.u^n , (1 + ^) II^W Ih 




(2.18) X ||A,CWIIi + IICWIIi + « 

V 

and 7 is a certain constant depending only on the properties of the 
domain Q. 

Observe that the dissipative solutions satisfy the initial condition 
(II. 6p . Indeed, at t = 0, inequality (12. ip becomes 

(2.19) 2fi\H0)-mfv+M0)-9{0)r < 2/i||a-C(0)||^+l|ao-^^(0)f 

This easily yields (II. 6p (see a similar reasoning below, in the proof 
of Proposition 14.11) . 

Moreover, these solutions obey the following "dissipative" estimate 

(2.20) 2^i\\u{t)\\l + \\a{t)f < 2fi\\u{0)\\l + ||a(0)f , Vt > 0. 

It follows from (11.61) and (12. ip with identically zero ( and 9. 

Remark 2.1. In the Euler- alpha case, the definition is much simpler: 
a function u G C^([0, 00); V) is called a dissipative solution to problem 
(O) - (fTTTTjl if, for all test functions C e C^{[0, 00); V3), and all non- 
negative moments of time t, one has 

\Ht) - ami < exp f / Tis)ds \\\\a- C(0)f^ 



MAXWELL-ALPHA FLUID 



9 



(2.21) +2 J exp TiOd^j (Ei(C, 0)(s), - C(s)) rfs} 
where 

r(t) = 7max{l, l/a'} (||A„CWI|i + IICWIli + "ICWIls) , 

and 7 is a certain constant depending only on Q. Obviously, (12.201) 
becomes 

(2.22) \\u{t)fy<\\uiO)fy, Vt>0. 

Our main result provides existence of dissipative solutions and their 
relation with the strong ones: 

Theorem 2.1. Let Q be a bounded domain having the cone property, 
a) Given a E V, ao €i L2, there is a dissipative solution to problem 

(HID - dH). 

b) If, for some a E V , ao & L2, there exist T > and a strong 
solution {ut,(Tt) e C\[0,T];V3) x C\[0,T];H^) to problem - 
(11.61) . then the restriction of any dissipative solution (with the same 
initial data) to [0,T] coincides with {ut,(Tt)- 

c) Every strong solution {u,a) G C"'^([0, 00); V3) x C^{[0, 00); H"^) is 
a (unique) dissipative solution. 

3. A REGULARIZATION AND PASSAGE TO THE LIMIT 

In order to prove Theorem 12.11 via approximation, we consider the 
following auxiliary problem: 

dv ^ dv ^ 
(3.1) -^ + ^^ Uig^ + ^Y1 + Vp + eAsu = 5 Div a, 

i=l * i=l 

^ + 5 J] Ui-^ + 6aW - 6Wa + eAicrj = 26r]£, 

(3.3) V = AaU, 

(3.4) divM = 0, 

(3.5) u =0, 

an 

(3.6) u\t=o = 5a, a\t=Q = 6ao. 
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Here, e > and < 5 < 1 are parameters. The first one will 
then go to zero, and the second one is needed for correct application 
of Schaeffer's theorem [121 P- 539]. We keep assuming a E V, ao E L2. 

The weak formulation of (13. ip - (13. 6 p is as follows. 

Definition 3.1. A pair of functions {u,a) from the class 

(3.7) ueWi, aeW2 

is a weak solution to problem (13. ip - (13. 6p if the equalities 



n 

i=l ' i=l 



(3.8) + £(n, <^)3 + 5(a, V<^) = 0, 
and 

-(.,*) + -(a,*) 

1=1 ^ 

(3.9) + 6{aW - Wa, $) + e{a, <l>)2 = 2(5/i(VM, $). 

are satisfied for all ^9 G V3, $ G //^(f2,Rg^") almost everywhere in 
(0,T), and (Q and hold. 

Remark 3.1. This notion of weak solution can be derived in a standard 
framework of multiplying by a test function, integrating by parts and 
using formula (12. 6p (see e.g. [H] concerning the Navier-Stokes system, 
see also [HI Section 6.1.1] for general remarks about weak formulation 
of equations). Note that formula (12. 5p should not be exploited at this 
stage. 

Lemma 3.1. Let Q be any domain having the cone property. Let 
u,a be a weak solution to problem (13.1 p - fl3.6p . Then, for all ( G 
C^{[0, 00); 1/3), e G C\[0, 00); H^) andO<t<T, one has 

2ii\\u{t) - mfv + \W{t) - 9{t)r 
t 

+2el{2f^Ms)-as)\\l + \\a{s)-e{s)\\l)ds 




<exp 5Tis)dsj |2/i||5a-C(0)ry+||5fTo-e(0)r 
exp ly 6T{0dn [Afi{E,{C,e,6){s),u{s)-as)) 
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+2(E2(C, 9, S)is), ais) - 9is)) - 4/i£(C(s), - C(s))3 

(3.10) -2eie{s),a{s) ~e{s))2]ds^ 
where T is as in (12.181) . 

Proof. Observe that 

J " n " 

-(C, - 5 5^(0A„C, ^) + 5 5^((A,C).V0, ^) + (i?i(C, 0, 5), if) 

i=l * 1=1 

(3.11) +e{C,ifh + 6{e,V<f) = e{C,^h, 
and 

|(^, $) + ji9, $) - 5 ^^(0^, ^) + 5(^w^(C) - W^(C)^^, $) 

i=l * 

(3.12) + {E^iC, e, 5), $) + £(^^, <l>)2 = 25/i(VC, $) + £(^, $)2. 

for 99 G Vs, $ G Denote w = u — C, and ^ = cr — For almost all 
t G (0,T), put if = w{t) and $ = ^(t). Multiply the difference between 
(13. 8 p and (13. lip by 2/i, and add this to the difference between (13. 9 p 
and (I3.12p . arriving at 

^d^^^'^^^+2d^^^'^) 

Tl r\ Tl r\ Tl r\ 

-2(3 /i 2_^{CAaW, — ) - 2^/i 2_^(wiA„C, g^) - 2a/i 2_^(w;iAc,w, — ) 

1=1 * i=l * i=l * 

71 n n 

+25/X ^((A,C)*Vty„ ti;)+25/i ^((A,u;)iVO, tf;)+25/i ^((A,ty),V«;i, w) 

i=l 1=1 1=1 

i=l i=l 

+S{9W{w) - Wiw)9, ^) + 5i<^Wiu) - W{u)<^, ^ + e{^, ^)2 + 2ne{w, w)^ 

(3.13) = 2/i(i?i(C, 9, 6),w) + {E^iC, 9, 5), ^ - 2fie{C, wh - e{9, 
We will need the following equalities: 

(3.14) - ^(KiA^K, ^) + ^(A„/€iV/€i, k) = 0. 

i=l * i=l 



(3.15) E('^^^'^) = 0- 



i=l 
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Here k G V3 and r G H^. The first equality is trivial. The second one 
is well-known and may be obtained via integration by parts. 
Note that 

n „ 

(3.16) -Wjk{t, x)<;ki{t, x)qji{t, x)) dx = 0, 
since ^ is a symmetric matrix. Moreover, 

(3.17) - Y,iw^^c.w, + Y^{{A^w),Vwi, w) = 0, 

i=l * i=l 

and 

(3.18) 5^(«.,, ^) = 



(by flXTi]) and fl315|) . respectively.). 

Now, let us estimate the remaining nonlinear terms in (13.13^ . Inte- 
grating by parts, we get 

i=l i=l *J=1 

W5^(C.v^„^). 

The first and the last terms vanish by (13.151) . The second one, due 
to the Cauchy-Bunyakovsky-Schwarz inequality and (12.111) . does not 
exceed Cna^CUMl By dZHD, 

^(u7iA„C, 7^) < Cn||A„C||i|kll?, 



dx 
1=1 



and 



Then, 



^{{/\^C)Nwi,w) < Cn\\A^C\\i\\w\\l 



i=l 



^((A„w;)iVCi,w) = ^(w,VO,w) 

i=l i=l 
ij=l i,j=l ^ 
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<Cn\\ai\\w\\l + Cna^(h\\w\\l 
Further, by integration by parts, 

1=1 

i=l i=l 

The first term is zero since w is divergence-free, and the second one, 
by f l232|) . is bounded by CnllwIlill^HalkH. Finally, by fl2lT]) . 

{6W{w)-W{w)6,c;)<Cn\\wU6h\\4- 
Now, for certain 7 depending on Q only, fl3.13p yields 



-(2i^\\w\\l + \\4^) + 2e{2fx\\w\\l+Ml] 



d 

di' 

< 25/iCn||A,C||i|kll? + 26fiCn\\C\\i\H\l 
+26fiCna'\\a4w\\l + 6Cn\\wU6hM 
+MEi{C 0, S),w) + 2{E2{C, 6, 5), - 4/i£(C, w)s - 2£(^, ^2 

< 57 (2/imax{l, l/an||u;||^(||A.Cl|i + llClli + "IClls + ll^lh) + lkl 

+MEi{C 0, S),w) + 2{E2{C, 6, 6),c;)- 4/i£(C, ^)3 - 2e{6, 
< 6T{2fi\\wfy + 

(3.19) + MEiiC, 0, S),w) + 2(E2(C, 0, 6),^)- 4/i£(C, wh - 2e{6, c;),, 

with r from (12.181) . It remains to apply Lemma 12.11 to this inequality. 

□ 

Lemma 3.2. Let u,a be a weak solution to problem (13. ip - (13.60 . The 
following estimates are valid: 

(3.20) 2^\\u\\l^(^o,T;V) + \W\\l^(o,T;L2) < 2/x||a||y + \\ao\\ = C, 

C 

(3-21) l|w|U2(0,T;y3) + \W\\l2{0,T;H'^) < 

(3.22) \\u\\L2iO,T;V*) + 1I^'IL,(0,T;H-) < ^^(1 + V^)- 

The constants C are independent of e and 5, but depend on \\a\\v, ||o"o||; 
T. 
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Proof. The estimates (13.201) and f l3.2ip are direct consequences of (13.11) 
with ( = 6 = 0. It remains to estimate the time derivatives, expressing 
them from (13.81) and (13. 9p . and taking into account (12. 5p and (12.60 : 



lUaCO.T) < S\\ ^{U^U, -^)\\l^(^o^t) + S\\ ^{UiVUi, ip)\\ L2{0,T) 
i=l * 1=1 



71 ^ n 

+6a^\\ ^{uiAu, ^)||l2(o,t) + 5a'^\\ ^{Aui\/Ui,(p)\\L2{o,T) 

i=l * 1=1 

(3.23) +e\\ {u, <^)3||l2(o,t) + S\\ (a, V(p) \\l2{o,t), 

and 

\\{(^',^)\\L2iO,T) < j\\{(^,^)\\l2{0,T) + S\\}_^{Uia, —)\\L,iO,T) 

1=1 * 



+6UaW-Wa,mL,io,T) 

(3.24) +£||(a,$)2|U2(o,T) + 2(5/i||(VM,$)lU2(o,T). 

Integration by parts implies 



dui du dip \ yr-^ ( du d'^Lp 

Ui 



^dxjdxj dxi J ^ V dxA dxidxj 
t,j=i ^ -J .J ' 2ij=i ^ -J -J 



(3.25) - ^ [ r^^\ , ] 

and 

n 



4 = 1 



dui dui d(pk\ / dui d'^u 



^ ^ \dxjdxk dxj) ' 9x,-9a;fc 

i,j,k=l ^ J ' i,j,k=l J J 

The second term is zero due to (I3.15p . 

Using Holder's and Sobolev imbedding inequalities, we estimate (jS]) 
as follows: 

W,ip)\\L2{Q,T) < ^||m||L(0,T;L2)I|V<^|U^ 
+^ll«IUoo{0,T;L2)||M||L2(0,T;y)||V(/?||L^ +25a2||^||2^^^^^^^jY(^||^^ 

+'^"l«IUoo(0,T;L4)l|w|U2(0,r;V)||VVl|L4 
+^l|w|U2{0,T;V3)IIV5||3 + '^l|cr||L2(0,T;L2)||Vv5||L2 
< C (||Ml|i^(o,T;y)ll^ll3 + ^l|M||L2(0,T;y3)ll^ll3 + 1^ lUoo (0,T;L2) II II s) 
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Now, (Km and ( ICTjl yield 

(3.27) ||(^/',<^)|U,(o,T)<C(l + v^)||v9||3, 

which is equivalent to the required bound for u'. 

Similarly, from (|3]) we derive the estimate for the time derivative of 

a: 

||(^',$)||l.(0,T) < C[\\a\\L,iO,T;L2)\M 
+ II^IUoo(0,T;L4)||cr||L2(0,T;L2)||V$||L4 + II Cr II L2(0,T;L2) || lUoc (0,r;l/) || * || 
+^I|o-||l2(0,T;H2)||'^'||2 + ll^^l|L2(0,r;y)||*||] 
< C'[lklUoo(0,T;L2)||'^'||2 + ||w||Loo(0,r;V)|k||Loo(0,T;L2)||'^'||2 
+^lklU2(0,T;H2)||$||2 + lh-f||Loo(0,T;y) ll'^'lh] 

(3.28) <C(l + v^)||$||2. 

□ 

Lemma 3.3. Given T > 0, a bounded domain Q having the cone 
property, and data a E V, ctq G L2, there exists a weak solution to 
problem f l3.ip - f l3.6p with 5 = 1. 

Proof. Observe that Lemma 13.21 yields the a priori estimate 

(3.29) \\u\\w, + \W\\w,<C, 

where C may depend on e but does not depend on 6. 

Let us rewrite the weak statement of f l3.ip - (13. 6p in the suitable 
operator form 

(3.30) A{u,a) = 5Q{u,a). 

For this purpose, we introduce the following operators. Here cp E V3 
and $ G are test functions. 

Ni:Wi^L2{0,T;V;), 

i=l 

n 
i=l 

Ns:WixW2^L2{0,T;H^^), 

$ 

dxi 



{Ns{u, a), $) = 2^(Mi(T, — ), 



1=1 
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N,:W^xW2^L2{0,T;H]^'), 
{N^{u, a), $) = {Wa - aW, $), 

S,:W2-^ L,{0,T;V;), {S^a,y^) = -(a, Vy^), 

S2:Wi^ L2(0, T; H^^), {S2U, $) = 2/i(Vn, $), 

53: 1^2^/^2(0, T; if/), '53a = -^, 

Q:W^xW2^ ^2(0, T; x ^2(0, T; ii/) x F x L2, 

Q(M,cr) 

= (Niiu) + N2iu) + Sia, Nsiu, a) + N^iu, a) + ^2^ + 5*3^, a, ctq), 

A:WixW2-^ ^2(0, T; V;) x ^2(0, T; if/) x x L2, 
A{u, a) = {u + eA^u, a' + £^2(7, u\t=Q, cr|t=o). 
We have treated the time derivative terms with the help of fl2.5l) and 

Let us check that the operator Q is compact (and, hence, continu- 
ous, since it consists of hnear, bihnear and quadratic components). It 
suffices to show compactness of 5*1, S'2, S'3, Ni, N2, N3 and A^4. 

Since the embeddings G V and C L2 are compact, the em- 
beddings Wi C 1^2(0, T; \/), W2 C 1^2(0, T; L2) are also compact by the 
Aubin-Simon theorem (see e.g. [39[|l8] or Theorem III. 2.1]). Thus, 
5*1, S'2 and S'3 are compact as superpositions of compact embeddings 
with bounded operators. 

Moreover, by [SHI Corollary 8], the embeddings 

C Lp{0,T;V),W2 C L,{0,T-L2), 

and 

W, C L^iO,T;V2) 

are compact for any p < 00 and g < 4. 

It remains to estabhsh boundedness (and, therefore, continuity) of 
the bilinear operators 

iV3 : L4(0, T; V) x L,{0, T; L2) ^ ^2(0, T; ii/). 



N, : L4(0, T; V) x ^4(0, T; L2) ^ ^2(0, T; fi'^). 
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N, : L,{0, T; V) x Ls{0, T; V2) ^ ^2(0, T; V,*), 



{Ni{u,w),(p) = ^{UiAaW, — ), 
i=l * 

N2 : ^^3(0, T; V^2) X ^6(0, T; V) ^ ^2(0, T; ^3*) 



i=l 



This is straightforward, e.g. for N2 one has 

II {N2{U, W), if) ||l2(0,T) < C\\A^u\\l;{0,T-M) II Vu^|Ue(0,T;L2) IIV^IUoo 
^ C'||M||L3(o,T;y2)II^IU6(0,T;y)||V^||3- 

The hnear operator A is invertible by [48i Lemma 3.1.3]. Thus, 
fl3.3ip can be rewritten as 

(3.31) {u,a) = 5A~^Q{u,a) 

in the space Wi x W2- 

By virtue of Schaeffer's theorem [121 p. 539], the a priori estimate 
fl3.29p . uniform in 5, guarantees existence of a fixed point of the map 
A^^Q, which is the required solution. □ 

We are now in a position to prove the main result. 

Proof. (Theorem 12.11) Take an increasing sequence of positive numbers 
Tra — 00 and a decreasing sequence of positive numbers — t- 0. By 
Lemma there is a pair (m^, a^) which is a weak solution to problem 
(13.11) - (13. 6p with (5 = 1, T = T^, e = Em- Denote by Um and am the 
functions which are equal to Um and am in [0, Tm] and are equal to zero 
on {Tm, +00). 

LemmaOimplies that, for all C G ^^([0, 00); 1^3), ^ e C\[0,ooy, H^) 
and < t < T < Tm, one has 

2fi\\um{t) - amv + \Wm{t) - dit)f 



<exp\^J T{s)dsj |2/i||a-C(0)rv^+||ao-^^(0)r 
+ jexp(j TiOdn [4f,{E,{C,e){s),Um{s)-C{s)) 
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+2(E2(C, 9)is), aUs) - Ois)) - Af^eUas),uUs) - C(s))3 

(3.32) -2€Ue{s), a^s) - ^(s))2] cis}. 

Fix an arbitrary interval [0,T]. Due to a priori estimate fl3.20p . 
without loss of generality (passing to a subsequence if necessary) one 
may assume that there exist limits u = lim Um, which is weak-* in 

Loo(0, oo; V) and weak in ^2(0, T; V), and a = lim a.^, which is weak- 

m—>-oo 

* in Loo(0, 00; L2) and weak in ^2(0, T; L2). 

Moreover, by (13.221) . without loss of generality one may assume that 
u'^ vl in L2(0,T; V^3*), a'^ a' in L2(0, T; i/^^). This gives that 
u e C{[0,T];V*), a e c([0,T]; H]^'^), and, by a well-known Lions- 
Magenes lemma, see e.g. [IBl Lemma 2.2.6], u G Cw{[0,T];V), a G 
a([0,T];L2). 

Take the scalar product in L2(0,T) of inequality ([3]) with a smooth 
scalar function ip with compact support in (0, T) and with non-negative 
values, and use (13.211) and the Cauchy-Bunyakovsky-Schwarz inequal- 
ity: 

T 

2^l\\Um{t) - C(t)||y + hm{t) - ^(t)ir}^(t) dt 


T 



< 





1 exp T{s)ds^ ^^2^^\\a-aml + \H-om'' 

+ y" exp l^y" r(Orfej [4/i(Ei(C,^^)(s),M^(s)-C(s)) 

(3.33) +2(E2(C, 0){s), a^{s) - e{s))] ds + Ci^ + em)}m dt. 

Passing to the limit inferior as m —t- 00 in ([3]), and using the fact 
that the norm of a weak limit of a sequence does not exceed the limit 
inferior of the norms, we arrive at: 

T 

I {2f^\\u{t) - cmv + Mt) - 9it)r]m dt 


T / t 

< jewij ns)ds I ^2/x||a-C(0)||^ + ||(7o-^^(0)|p 
\o 
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+2(E2(C, 9){s), ais) - e{s))] ds^i^it) dt. 




(3.34) 



Since ip and T were chosen arbitrarily, ([3]) yields fl2.ip . and we have 
proven the existence of a dissipative solution. Note that (12. ip holds at 
every non-negative moment of time owing to the weak continuity of u 
and a. 

Now, let {u, a) be a dissipative solution with the same initial data 
as the strong solution {uT,crT)- Putting ( = ut, 9 = ar in (12. ip for 
t e [0,T], and taking into account that Ei^ut.o^t) = E2{uT,crT) = 
on [0,T], we get that the right-hand side of (12. ip vanishes there, and 
we arrive at the claim b) of Theorem 12.11 Observe that c) is a direct 
consequence of a) and b): any strong solution should coincide with all 
dissipative solutions as long as it exists (cf. [45j). □ 

Remark 3.2. Remember that, in Lemmas 13.11 and 13.21 Q can be an 
unbounded domain with the cone property. Moreover, the constants 
in (I3.20p and (I3.2ip do not depend on Q. Then, in order to establish 
existence of a dissipative solution to (II. ip - (II. 6p . one can try to ap- 
proximate Q with regular bounded domains (cf. e.g. [43j) and pass to 
the limit, without generalizing Lemma [3.31 onto the case of unbounded 
ones. However, some pitfalls appear, in particular, we did not manage 
to prove the weak continuity of the limiting function. So we leave the 
general "unbounded" case as an open problem. In the particular case 
fl = W^, Lemma [3.31 and Theorem 12.11 are valid since one can replace 
the operators and A2 and the spaces V3 and with (/ — A)^, 
(/ — A)^, the closure of V in H^, and Hq, respectively, and approxi- 
mate Q with an ascending sequence of concentric balls (see a similar 
reasoning in [15]). 

Remark 3.3. Assume again a G V, ctq € I/2. Take a sequence of 
positive numbers am — ^ 0. The dissipative bound (I2.20p guarantees 
compactness of any sequence {{um, crm)} of dissipative solutions to the 
Maxwell-am, problem in the weak-* topology of Loo{0, oo; H x L2). 
The accumulation points of {{umyCTm)} can be considered as "ultra- 
generalized" solutions to the IBVP for motion of the Maxwell fluid, i.e. 
the problem (II. ip - (I1.6P with a = 0. However, the relevance of this 
notion is an open question. 



Remark 3.4. A related open problem is whether our dissipative solu- 
tions to the Euler-a problem converge in some sense to the dissipative 
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solutions to the Euler problem as a — )■ 0. The convergence is known 
for the strong solutions on the time interval where they exist 

Remark 3.5. Wild solutions [lOj to the Euler equations are an im- 
portant example of dissipative solutions. These solutions are known 
to exist only for rough initial data and to be non-unique. A slightly 
different sort of wild solutions, which are not dissipative solutions, but 
exist for all initial data, was considered in [IS]. It is interesting whether 
the a-models admit wild solutions. 

4. Appendix. Dissipative solutions to Cauchy problems in 

HiLBERT SPACES 

In this appendix we illustrate the idea of dissipative solution by con- 
sidering (a little bit informally) this notion for an abstract differential 
equation in a Hilbert space. We point out that this is not a universal 
definition of dissipative solution but just a guideline, which shows the 
essence of this concept. Slight natural modifications of this approach (if 
necessary) may lead to definitions of dissipative solutions for particular 
PDEs, e.g. the corresponding definitions from |26lB5] or Definition 12.11 
of the current paper. 

Let X be a Hilbert space with inner product (■, ■) and Euclidean 
norm || ■ ||, and F : [0, -foo) x X — )■ X be a nonlinear operator, which 
can be discontinuous. Assume that F satisfies the following condition: 

(4.1) (Fit, x) - Fit, y),x-y)< d{t, y) \\x - , Vx, y G X, t > 0, 
with some locally bounded function d{t, y) of t > and ?/ G X. 
Remark 4.1. For example, (14.1 p holds provided 

(4.2) {F{t, x), x) < 0, Vx G X, t> 0, 

and the operator F can be decomposed into a sum of a linear (in x) 
operator A and a quadratic (with respect to x for fixed t) operator 
subject to a slight boundedness assumption. Indeed, let 

F{t,x) = A{t,x) + /(t,x,x) 

be such a decomposition, where / is a bilinear operator with respect 
to (x, y) satisfying the bound 

||/(t,x,y)|| <c(t)||x||l|i/|| 

with a locally bounded function c. Take any x, y G X, and let z = x—y. 
Then 

(F(t, x) - F{t, y),z) = {F{t, z),z) + (/(t, y, z),z) + (/(t, z, y), z) 

<2c{t)\\ym?- 



MAXWELL-ALPHA FLUID 



21 



We consider the abstract Cauchy problem 

(4.3) u'{t) = F{t,u{t)), t > 0, 

(4.4) m(0) = aeX. 

Let TZ be some fixed set of differentiable functions from [0, +oo) to 
X. The solutions u E 71 to (14. 31) . (14. 4p will be called (sufficiently) 
regular. More generally, a solution u : [0,T] — )■ X, T > 0, is said to be 
regular if it coincides with the restriction of a function from TZ to this 
interval [0, T]. We assume that 71 is sufficiently large, so that the set 
{u{0)\u G 7Z} is dense in X. 

Take a solution u e7Z and any test function v eTZ, and denote 

E{t,v{t)) = -v\t) + F{t,v{t)). 

Then (14. 3 p implies 

(4.5) {u - v)' = F{t, u) - F{t, v) + E{t, v). 
Hence, 

[\\u-vfy = 2{u'-v',u-v) = 2{F{t,u)-F{t,v),u-v)+2{E{t,v),u-v) 

(4.6) < 2d{t,v{t))\\u - + 2{E{t,v),u - v). 
An application of Lemma 12.11 yields 

u{t) - v{t) f < exp I J 2d{s, v{s))ds 

.0 



2d{^, viO)d^ I (Eis, vis)), uis) - v{s)) ds 

Definition 4.1. A weakly continuous function u : [0, +oo) — )■ X is 
called a dissipative solution to problem (l4.3p .f H^ provided (j4]) holds 
for all f G 7^ and t > 0. 

The basic properties of these solutions are given by 

Proposition 4.1. a) If, for some a E X, there exist T > and a regu- 
lar solution ut '■ [0,T] X to problem (l4.3p .f H^ . then the restriction 
of any dissipative solution (with the same initial data) to [0,T] coin- 
cides with Ut- b) Every regular solution u e7Z is a (unique) dissipative 
solution, c) Any dissipative solution satisfies the initial condition (14. 4p . 
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Proof. Let -u be a dissipative solution with the same initial data as Mr- 
There exists ut ^Ti coinciding with ut on [0, T]. Putting f = in (jl]) 
for t G [0, T], and taking into account that E{t, Urit)) = on [0, T], we 
conclude that the right-hand side of (jlj) vanishes there. Thus we have 
a). The reasoning above which led us to the derivation of inequality 
(jlj) gives us b); note that the uniqueness in b) follows from a). To show 



c), we put t = in (jl]) and get 

(4.8) ||m(0) -t;(0)|| < ||a-t;(0)||. 
Since the set of possible f (0) is dense in X, we have 

(4.9) ||m(0)-6|| < ||a-6|| V6gX, 



and it remains to let 6 = a. □ 

Remark 4.2. The uniqueness of dissipative solutions is not given a 
priori but follows from the existence of a regular solution. On the 
other hand, when a regular solution does not exist, non-uniqueness is 
possible (see Remark [3.51) . 

Remark 4.3. Consider the "quasihomogeneous" case, i.e. when 

(4.10) F(-,0)=0. 
If 

(4.11) rf(-,0) = 0, 

then (|4]) with v = Q yields dissipative behaviour of solutions, 

(4.12) ||M(t)|| < ||m(0)||, t > 0. 

All this holds true, in particular, in the framework of Remark 14. 1[ A 
similar situation happens for the homogeneous 3D Euler equation [3]. 
However, Definit ion 14 . 1 1 may be useful without assumptions of this kind 

(cf. m)- 

A remarkable thing is that these dissipative solutions exist under 
minimal assumptions on F. We are not going to formulate an existence 
theorem, but we'll try to explain why they exist. 

The first point is an a priori estimate. The Cauchy-Bunyakovsky- 
Schwarz inequality and (jH) with v = yield 

(4.13) {\\ury<2(^d{t,0) + \^ \\ur + 2\\F{tM\\ 
so, by Lemma 12.11 we have the following bound: 
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(4.14) X \\a\\' + 2 exp\ 2[d{i,0) + -]dn\\F{s,0)\\'ds 



Now approximate F by smooth functions F^, e > 0, F^{-,v{-)) 
F(-,t>(-)) in Li ;oc(0, oo) for any fixed f G 7?. as e — 0. The functions 
F^ satisfy (14. ip with some (since they are smooth), and thus the so- 
lutions of the corresponding systems fl4.3p .( H^ are a priori bounded 
by analogues of ([1]). The infinite-dimensional Picard-Lindelof theorem 
(see e.g. [321 137]) implies that the solutions exist (globally) and are 
unique. Assume that there exists a function d{t, y) for which (14. ip is 
valid for any e > 0. Without loss of generality, d = d {if not, one 
can replace both of them by max[c/, ). If the functions F^{-,0) are 
locally square integrable in time, and the corresponding integrals are 
uniformly (i.e. independently of e) bounded, then the uniform a priori 
estimate holds for u^. Thus, any sequence {u^,.}, e/c 0, is rela- 
tively compact in the weak-* topology of -Loo(0, T; X) for any T > 0, so 
without loss of generality there exists a limit u. A diagonal argument 
guarantees that the sequence and the limit can be chosen so that they 
do not depend on T. 

Suppose that there "large" Banach space Y (so that X is 

continuously embedded into it), and, say, a continuous function c : 
^ M such that 

\\F,{t,x)\\Y < c{t, \\x\\x), xeX, t>0. 

This assumption implies a uniform bound for the time derivatives 

lkelUoo{0,T;y) < C. 

Thus, u' G Loo(0,T;F), and, by HH Lemma 2.2.6], u G C^([0,T];X). 

Now, remember that u^^, are, in particular, dissipative solutions. 
Multiplying by a smooth nonnegative function tp{t) and integrating 
from to T, we can pass to the limit in (jl]) in the same way as it 
was done in the proof of Theorem 12. 1^ and we conclude that m is a 
dissipative solution to the original system f l4.3p . fl4.4p . 

Remark 4.4. Let X be finite-dimensional. Then one can study system 
(|0]),(|lill) via Filippov's approach [El E]. Condition (gl]) locally 
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implies the so-called right Lipschitz condition 

{F{t,x)-F{t,y),x-y)<C\\x-y\\\ 

which guarantees uniqueness of Filippov solutions. On the other hand, 
it is not clear whether the Filippov solution is always a dissipative one, 
or vice versa (or at least whether dissipative solutions are unique in 
this case). 

Remark 4.5. The model considered recently in [40| does not satisfy 
any condition similar to (14. ip . but we found an analogue of (jl]). How- 
ever, that inequality has an absolute value in the last integral, so the 
weak passage to the limit is not an option. We proceed there via some 
strong convergence, which is still not strong enough to get classical 
weak solutions. 
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